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CN . ABSTRACT 



, Context. The majority of chemical evolution models assume that the Galactic disk forms by means of infall of gas and divide the disk into 
1^ ■ several independent rings without exchange of matter between them. However, if gas infall is important, radial gas flows should be taken into 
\^ account as a dynamical consequence of infall. 
^ Aims. The aim of this paper is to test the effect of radial gas flows on detailed chemical evolution models (one-infall and two-infall) for the 
^' Milky Way disk with different prescriptions for the infall law and star formation rate. 
I Methods. We modified the equation of chemical evolution to include radial gas flows according to the method described in Portinari & Chiosi 
O Results. We found, that with a gas radial inflow of constant speed the metallicity gradient tends to steepen. Taking into account a constant time 
^ ' scale for the infall rate along the Galaxy disk and radial flows with a constant speed, we obtained a too flat gradient, at variance with data, 
^ , implying that an inside-out formation and/or a variable gas flow speed are required. To reproduce the observed gradients the gas flow should 
■ ' ■ increase in modulus with the galactocentric distance, both in the one-infall and two-infall models. However, the inside-out disk formation 
coupled with a threshold in the gas density (only in the two-infall model) for star formation and/or a variable efficiency of star formation with 
galactocentric distance can also reproduce the observed gradients without radial flows. 
^ ■ Conclusions. We showed that the radial flows can be the most important process in reproducing abundance gradients but only with a variable 
gas speed. Finally, one should consider that uncertainties in the data concerning gradients prevent us to draw firm conclusions. Future more 
detailed data will help to ascertain whether the radial flows are a necessary ingredient in the formation and evolution of the Galactic disk and 
disks in general. 
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1, Introduction are in very good agreement with the data on Cepheids in the 

galactocentric distance range 5-17 kpc for many of the studied 



■ Generally, a good agreement between model predictions and eje^ients 



observed properties of the Galaxy is obtained by models which 
c3 assume that the disk formed by infall of gas: (e.g Matteucci & Colavitti et al. (2009) found that it is impossible to fit all the 
Francois (1989), Chiappini et al. 1997, Frangois et al. 2004). disk constramts at the same time without assuming an inside- 
The formation timescale of the thin disk is assumed to be a formation for the Galactic disk together with a threshold 
function of the galactocentric distance, leading to an inside-out *e gas density for the SF rate. In particular, the inside-out 
scenario for the Galaxy disk build-up (Matteucci & Frangois formation is important to reproduce the right slope of the abun- 
1989). This assumption helps in reproducing the abundance ^ance gradients in the inner disk, whereas the threshold gas 
gradients along the disk as shown in studies such as Chiappini density is important to reproduce present day gradients in the 
et al. (2001), Cescutti et al. (2007), Colavitti et al. (2009) and o^te^" ^^i^^. Generally, simple models without radial gas flows 
Marcon-Uchida et al. (2010). In these papers were analyzed ^ ""^dially constant formation timescale cannot repro- 
the parameter space of the "static" chemical evolution models ^uce the slope in the inner disk. In the framework of models 
considering diff-erent infall and star formation (SF) rate laws in no inside-out mechanism such as the model with the cos- 
order to reproduce the observed abundance gradients. Cescutti ^ological infall law of Colavitti et al. (2009) they tested the 
et al. (2007) showed that the results obtained using a chemical ^^^""^ *e efficiency of SF varying with the galactocentric 
evolution model with a two-infall law and inside-out formation ""^dius, being higher in the innermost than in the outermost 

regions of the Galactic disk. This assumption, even without 

* email to: spitoni@oats.inaf.it the threshold in the surface gas density, can produce gradients 
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with the right slope both in the inner and outer regions of the 
Galactic disk but it fails in reproducing the gas density distri- 
bution along the disk. However, if gas infall is important, the 
Galactic disk is not adequately described by a simple multi- 
zone model with separated zones (Mayor & Vigroux 1981). To 
maintain consistency, radial gas flows have to be taken into ac- 
count as a dynamical consequence of infall. The inf ailing gas 
has a lower angular momentum than the circular motions in the 
disc, and mixing with the gas in the disc induces a net radial 
inflow. Lacey & Fall (1985) estimated that the gas inflow ve- 
locity is up to a few km s"^ Goetz & Koppen (1992) studied 
numerical and analytical models including radial flows. They 
concluded that radial flows alone cannot explain the abundance 
gradients but are an efficient process to amplify the existing 
ones. In particular, they suggested that the speed of the flow 
should increase with galactocentric distance and only change 
by 0.15 km s"^ per kpc, well below any detection limit for drift 
velocities of molecular clouds. Portinari & Chiosi (2000) im- 
plemented in the detailed chemical evolution model of Portinari 
& Chiosi (1999), characterized by a single infall episode, the 
eff'ects of radial gas inflows by means of the following assump- 
tions: 

- an uniform inflow pattern through the disk, with speed yr = 
-1 km s"^ applied to a model with a Schmidt (1959) SF law 
and they found , as expected, that the metallicity gradient 
tends to steepen in the presence of radial gas flows; 

- they considered models with various SF laws and tuned the 
inflow velocity pattern, in each case, as to match the obser- 
vational data relative to the radial profiles of the Galactic 
disk. They suggested that inclusion of radial flows in the 
chemical models can improve the match with the data. For 
these models the radial velocities span the range -1 < yr 
< -0.3 km s-^ 

In the work of Schonrich & Binney (2009) both the stellar and 
gas flow were considered in a chemical evolution model with- 
out inside-out formation. They are able to reproduce the ob- 
served data considering a variable speed of the gas flow in the 
range km s"^ to -5 km s"^ 

Another physical cause of radial flows is the gas viscos- 
ity. Viscosity in the gas layers induces radial inflows in the 
inner parts of the disc and outflows in the outer parts, with 
velocities of -0.1 km s"^ (Lacey & Faff 1985). Thon & 
Meusinger (1998) presented a model for the chemical and dy- 
namical evolution of the Galactic disk which combines viscous 
radial gas flows with infall of external gas onto the disk, and 
inf all-induced radial gas flows. In this hydrodynamic "hybrid" 
model it was assumed that the viscosity of the gas is closely 
related to the processes of SF and the energetic response from 
massive stars, i.e. the timescales for viscous transport and for 
gas consumption by SF are nearly the same. Following the pa- 
rameterization of Sommer-Larsen & Yoshii (1989), they ob- 
tained a SF rate law depending on the viscosity parameters. 

Gravitational interactions between gas and spiral density 
waves lead to large-scale shocks, dissipation and therefore ra- 
dial inflows of gas (or outflows in the outer parts) with typical 
velocities of -0.3 km s"^ (e.g. Berlin & Lin 1996 and refer- 
ences therein); much larger velocities can be achieved in the 



inner few kpc in the presence of a barred potential. However, 
given the small values of the flow speed in these works, the 
efl'ect on the abundance gradients are negligible. 

In this work we include radial inflows in 2 difl'erent chem- 
ical evolution models: a one-infall and a two-infall model. In 
particular, this is the first time that radial flows are included in 
the two-infall model. We test a constant velocity pattern for the 
radial inflow as well as a best fit pattern as done in Portinari & 
Chiosi (2000). We then study the eff'ect of radial flows on the 
formation of gradients coupled with the inside-out mechanism 
formation and the eff'ect of the threshold in the SF. Last but not 
least, we compare our results with the most recent data now 
available, including Cepheids,and this represents a novelty rel- 
ative to previous works. This comparison will allow us to put 
constraints on the mechanisms of gradient formation. 

The paper is organized as follow: in Sect. 2 we describe the 
reference "static" models used in this work, in Sect. 3 we re- 
port the nucleosynthesis prescriptions, in Sect. 4 we present the 
implementation of the radial inflow of gas in a detailed chemi- 
cal evolution model, in Sect. 5 observational data are shown. In 
Sect. 6 we report and discuss our results. Finally we draw the 
main conclusions in Sect. 7. 

2. The chemical evolution model for the Milky Way 

In order to reproduce the chemical evolution of the thin-disk, 
we adopt as a reference model an updated one-infall version 
of the chemical evolution model presented by Matteucci & 
Fran9ois (1989). In this model, the galactic disk is divided into 
several concentric rings which evolve independently without 
exchange of matter. 

The infall law for the thin-disk is defined as: 

B(rj) = b(r)e~^, (1) 

where td is the timescale for the infalling gas into the thin- 
disk. It is worth noting that eq.O diff'ers from the two-infall 
model of Chiappini et al. (1997) which assumes an infall law 
such as: 

A(r, = a{r)e-"'^^'^ + j^^^^^-it-t^a.yroir) ^ ^2) 

To have an inside-out formation on the disk, the timescale 
for the mass accretion is assumed to increase with the Galactic 
radius following a linear relation given by (see Chiappini et al. 
1997): 

Td = 1.033r(kpc) - 1.27 Gyr, (3) 

for galactocentric distances > 4 kpc. The region within 4 kpc 
is not considered in this paper and has by the strict inflow no 
impact on the outer regions. The coefficient b(r) is obtained by 
imposing a fit to the observed current total surface mass den- 
sity in the thin disk as a function of the galactocentric distance 
given by: 

cr(r) = CToe-'''^ (4) 

where (T/)=531 M© pc"^ is the central total surface mass den- 
sity and r/) = 3.5 kpc is the scale length. 
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In order to make the program as simple and generalized as 
possible, we used a SF rate proportional to a Schmidt (1959) 
law: 



jA(r,Oocyi:J^,(r,0 



(5) 



where v is the efficiency in the SF process and the surface 
gas density is represented by ^gas(r, t) while the exponent k is 
equal to 1.4 (see Kennicutt 1998). 

The equation below describes the time evolution of G/, 
namely the mass fraction of the element / in the gas: 

Gi{rj) = -ilf{rj)Xi{rj) 



MBm 

+ J" i//(r,t-Tm)Qmi(t-Tm)(f>(rn)dm 

Ml 

Mbm 
Mb, 



Mbm 0.5 

J (P(Mb) • J J 



dMn 



Mbm 

+ (1 - Ala) J t - Tm)Qmi{t - Trn)(p(m)dm 

MBm 



Mn 



f J" il/{rj- 

Mbm 



Tm)Qmiit - Tm)(p(m)dm 



(6) 



where Gi(r,t) = [crg(r,t)Xi(r,t)]/crA(r), crg(r,t) is the surface 
gas density, and (TA(r) is the present- time total surface mass 
density. Xa- are the abundances in the infalling material, which 
is generally assumed to be primordial. X/(r, t) is the abundance 
by mass of the element / and Qmi indicates the fraction of mass 
restored by a star of mass m in form of the element /, the so- 
called "production matrix" as originally defined by Talbot & 
Arnett (1973). We indicate with Ml the lightest mass which 
contributes to the chemical enrichment and it is set at O.SM©; 
the upper mass limit, Mu, is set at lOOM©. 

For the IMF, we use that of Scalo (1986), constant in time 
and space. Tm is the evolutionary lifetime of stars as a function 
of their mass m (Maeder & Meynet 1989). 

The Type la SN rate has been computed following Greggio 
& Renzini (1983) and Matteucci & Greggio (1986) and it is 
expressed as: 



RsNela 



Mbm r 0.5 

■ Ala J (P(Mb) J J 



Mb, 



f(jl)l//(t -TM2)djd 



dMn 



(7) 



where M2 is the mass of the secondary, Mb is the to- 
tal mass of the binary system, ji = M2/MB, iim = 
max[M2(t)/MBAMB-0.5MBM)/MBl Mbm = ^Mq, Mbm = 
I6M0. The IMF is represented by (P(Mb) and refers to the total 



mass of the binary system for the computation of the Type la 
SN rate, f(jj,) is the distribution function for the mass fraction 
of the secondary: 



/0u) = 2i^^(l+r)yu^ 



(8) 



with y = 2; Ala is the fraction of systems in the appropriate 
mass range, which can give rise to Type la SN events. This 
quantity is fixed to 0.05 by reproducing the observed Type la 
SN rate at the present time (Mannucci et al. 2005). Note that in 
the case of the Type la SNe the "production matrix" is indicated 
with 2^f^^ because of its diff'erent nucleosynthesis contribu- 
tion (for details see Matteucci & Greggio 1986 and Matteucci 
2001). 

We also apply the eff'ect of radial inflows of gas on the 
model described in Francois et al. (2004). 

We do not enter into details of this model, but we just 
remind here that in this case the Galaxy is assumed to have 
formed by means of two main infall episodes: the first formed 
the halo and the thick disk, the second the thin disk. (i.e. two- 
infall model) The typical time- scale for the formation of the 
halo is 0.8 Gyr and the entire formation period for the halo- 
thick disk does not last more than 2 Gyr. The time-scale for the 
thin disk is much longer, 7 Gyr in the solar vicinity, implying 
that the infalling gas forming the thin disk comes mainly from 
the intergalactic medium and not only from the halo (Chiappini 
et al. 1997). Moreover, the formation timescale of the thin disk 
is assumed to be a function of the galactocentric distance, lead- 
ing to an inside-out scenario for the Galaxy disk build-up (eq. 
O. The galactic thin disk is approximated by several indepen- 
dent rings, 2 kpc wide, without exchange of matter between 
them. 

The main characteristic of the two-infall model is an almost 
independent evolution between the halo and the thin disk (see 
also Pagel & Tautvaisienne 1995). A threshold gas density of 
IMqPc'^ in the SF process (Kennicutt 1989, 1998, Martin & 
Kennicutt 2001) is also adopted for the disk. We consider for 
the halo a constant surface mass density as a function of the 
galactocentric distance at the present time equal to 17 M© pc~^ 
and a threshold for the halo phase of 4 M© pc~^ as it was as- 
sumed for the Model B in Chiappini et al. (2001). 

The time-steps are small for the first 2 Gyr, of the order 
of 10^ - 10^ years, then they become larger according to the 
smaller variation of the variables and reach values of the order 
of 1 Gyr. 

3. Nucleosynthesis prescriptions 

For the nucleosynthesis prescriptions of the O and Fe we adopt 
those suggested in Fran9ois et al. (2004). They compared theo- 
retical predictions for the [el/Fe] vs. [Fe/H] trends in the solar 
neighborhood for the above mentioned elements and they se- 
lected the best sets of yields required to best fit the data. In 
particular they found that the Woosley & Weaver (1995) metal- 
licity dependent yields of SNe II provide the best fit to the data. 
No modifications are required for the yields of Ca, Fe, Zn and 
Ni as computed for solar chemical composition. For oxygen, 
the best results are given by the Woosley & Weaver (1995) 
yields computed as functions of the metallicity. For the other 
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elements, variations in the predicted yields are required to best 
fit the data (see Francois et al. 2004 for details). 

Concerning O, the best agreement with the [0/Fe] vs. 
[Fe/H] and the solar O abundance, as measured by Asplund 
et al. (2009), is obtained by adopting the original Woosley & 
Weaver (1995) yields from massive stars as functions of metal- 
licity. The same is not true for Fe: it was already known from 
the paper of Timmes et al. (1995) that the Fe yields as func- 
tions of metallicity overestimate the solar Fe and many people 
use those yields divided by a factor of 2. Alternatively, one can 
use the yields for solar chemical composition for the whole 
galactic lifetime and the result is the same, and this is what we 
do here. 

Concerning the yields from Type SNela, revisions in the 
theoretical yields by Iwamoto et al.(1999) have been taken into 
account.The prescriptions for single low-intermediate mass 
stars are by van den Hoek & Groenewegen (1997), for the case 
of the mass loss parameter, which varies with metallicity (see 
Chiappini et al. 2003, model5). 

4. The implementation of the radial inflow 

Lacey & Fall (1985) showed that several mechanisms could 
drive a radial gas flow in the Galactic disk. We refer here to 
the possibility that the specific angular momentum of the in- 
falling gas diff'ers from the value for the circular motions in the 
disk. Therefore, mixing will induce a radial flow. Lacey & Fall 
(1985) suggested that a good estimate for velocity of that gas at 
10 kpc is V/? = - 1 km s"^ We implement radial inflows of gas 
on our reference model following the prescriptions described 
in Portinari & Chiosi (2000). 

Let the ^-th shell be defined by the galactocentric radius 
rk, its inner and outer edge being labeled as r^_i and r^+i. 
Through these edges, gas inflow with velocity v^_i and v^+i, 
respectively. The flow velocities are taken positive outward and 
negative inward. 

Radial inflows, with a flux F(r), contribute to alter the gas 
surface density in the ^-th shell in the following way: 



where jSk and yk are, respectively: 



dcTn 



dt 



1 



rf 



[F(r,,i)-F(r,_pl 



where the gas flow at 1 can be written as: 



(9) 



(10) 



As in Portinari & Chiosi (2000) we take the inner edge 
of the shell, r^_i, at the midpoint between the characteris- 
tic radii of the shells k and k - I, and similarly for the outer 

edger^^i: 

r^_i = (rk-i + rk)/2, and r^+i = (rk + a+i)/2. We get: 



I 2 2 \ n+i- rk-i I 



rk-i + rk+i 



Inserting eqs. ([TOl) and (fTTI) into eq. ©, we obtain the radial 
flow term to be added in eq. ©: 



d 
dt 



(12) 



'k ^ 2 



rk-i + rk 
2 rk+i - Vk-i 



Jk 



rk + rk+i 



A{k+l) 
O-Ak 



(13) 



(14) 



rf 



crA{k+\) and CTAk are the actual density profile at the radius 
Vk+i and Vk, respectively. This formulation for the radial inflow 
follows the one used in Portinari & Chiosi (2000) but we want 
to stress that in the quantity Pk, concerning the term in square 
brackets, contains Vk+i - Vk-i instead of Vk+i + Vk-i as reported in 
the work of Portinari & Chiosi (2000), but certainly they used 
the right expression in their calculations. We assume that there 
are no flows from the outer parts of the disk where there is no 
SF. In our implementation of the radial inflow of gas, only the 
gas which resides inside the Galactic disk within the radius of 
18 kpc can move inward by radial inflow. Our model is not 
subject to numerical instabilities up to 18 kpc. Therefore since 
we choose not to consider flows from galactocentric distances 
larger than 18 kpc, our results are not dependent on these insta- 
bilities. 

5. Observational data 

For the galactic abundance gradient of oxygen we use the fol- 
lowing set of data: Deharveng et al. (2000), Esteban et al. 
(2005), Rudolph et al. (2006), who analyzed the Galactic HII 
Regions; Costa et al. (2004) who studied Planetary Nebulae 
(PNe); Andrievsky et al (2002a,b) who analyzed Galactic 
Cepheids. In order to compare our results with homogeneous 
sets of data, we separate the data by HII regions and PNe from 
the Cepheid ones. To better understand the trend of the data, 
we divide the data in 6 bins as functions of the galactocen- 
tric distance. In each bin we compute the mean value and the 
standard deviations for the studied elements. In Fig. [T] we re- 
port the whole collection of the data used in this paper. We 
are aware that there are some severe systematic uncertainties in 
each study that might justify larger error bars. 

Concerning the data for the surface gas density profile and 
the SF rate along the Galactic disk we use the Dame et al. 
(1992) and Rana et al. (1991) ones. 

6. The model results 

In this section the results of both the two-infall model of 
Chiappini et al. (1997) and the one-infall model developed in 
this paper are presented. The two infall model can have an ef- 
fect on the disk evolution but only at large galactocentric dis- 
tances, as discussed in Chiappini et al. (2001). 

We consider that in all models the efficiency of the SF is 
fixed at the value of v = 1 Gyr"^ in agreement with previous 
papers and in order to minimize the variation of model param- 
eters. 

In this paper we will focus on the eff'ect of radial flows with 
difl'erent prescriptions for the infall and the SF rate. In Table [T] 
we summarize all the models we consider as a function of the 
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Fig. 1. Radial abundance gradient for oxygen from observa- 
tions. The data are taken by Deharveng et al. (2000) (filled 
blue squares), Rudolph et al. (2006) (filled green pentagon), 
Costa et al. (2004) (red filled triangles), Andrievsky et al. (2002 
a,b) (magenta filled hexagons), and Esteban et al. (2005) (black 
filled circles). We label with the large blue squared points the 
mean Cepheids values and relatives error bars, whereas with 
the large black filled circles the mean values form the HII and 
Planetary Nebulae regions. 
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Fig. 2. Radial abundance gradient for oxygen. The black solid 
line refers to the one-infall model without threshold (SI I), the 
blue long dashed line to a one-infall model with a threshold of 4 
MqPc~^ (S1IT4) and the red dotted line a model with a thresh- 
old of 7 MqPc~^ (S1T7). The filled circles and relatives error 
bars are the observed values from HII regions and Planetary 
Nabulae. 

type of infall, the timescale of the thin disk phase, the presence 
of the threshold, and the type of radial inflow. In particular, 
in the first column we report the model names, in the second 
the infall type, in the third in the forth it is specified if there 
is inside-out formation (td variable with the radius) and the 
threshold in the SF, respectively. Einally, in the fifth it is in- 
dicated the type of the implemented radial flow (only for the 
"no- static" models). 




Fig. 3. Radial abundance gradient for oxygen. The black solid 
line refers to the one-infall model without threshold and with 
a constant = 4 Gyr (SIN). The blue long dashed line rep- 
resents the model with a radial inflow of 1 km/s (RlNl). The 
data are the same of Eig. [2l 



Eirst of all, we present the results obtained with the one- 
infall model with an inside-out formation. In Eig.|2]we show 3 
difl'erent cases without radial flow: a model without threshold 
in the SE (SH), models SnT4 and S1T7 with a gas threshold 
for SE of 4 and 7 M© pc~^, respectively, compared with the 
data from HII and PNe. We conclude that even if we consider a 
model with an inside-out formation the obtained gradient with 
a one-infall model without threshold is too flat and the obser- 
vational data are not reproduced especially in the outer part of 
the Galaxy disk. We see instead that the slope for galactocen- 
tric distances between 4 and 14 kpc can be reproduced if we 
consider the model with a threshold in the SE. However, since 
the surface mass gas density in the outer parts of the disk is too 
small, there is no SE in this region and then no metals produc- 
tion and no chemical evolution. 

Therefore, a possible solution to reproduce the data can be 
to consider a radial inflow of gas. 

Portinari & Chiosi (2000) computed the efl'ects of radial 
flows on their one-infall model without an inside-out forma- 
tion. They adopted an unique timescale for disk formation irre- 
spective of the galactocentric distance. Surprisingly, they found 
also in the case without radial flows an abundance gradient with 
a significant slope along the disk for the oxygen. 

In the Eig. [3] we report the model SIN using the same pre- 
scriptions of the SIl model but without the inside-out forma- 
tion, namely keeping td constant along the Galaxy disk and 
equal to 4 Gyr. We find a flat gradient, and in the outer part 
of the disk the values of log(0/H)-hl2 are even larger than the 
value in the inner part of the Galaxy. Portinari & Chiosi (2000) 
labeled SI 5a a model with the follow main prescriptions: i) 
one-infall, ii) no threshold, iii) no inside-out. With these pre- 
scriptions they found a significant gradient along the Galactic 
disk of ~ -0.03 dex kpc"^ in contrast with our model result 
reported in Eig. [3] with the black solid line. 

In Eig. [3] we also show the model RlNl results where it was 
considered a constant speed of 1 km s"^ for the radial inflow on 
the one-infall model no threshold and no inside-out formation. 

We see that the main efl'ect of a migration of gas from the 
outer part of the Galaxy toward the inner part without thresh- 
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Table 1. The list of the models described in this work. The first letter in the model names is S if that model is a "static" one, 
otherwise R if we show a model with radial flow of gas. The follow number (1 or 2) is related to the number of infall, the letter I 
is for the case for inside-out model and N for a constant formation time scale along the disk. When we consider a threshold in the 
star formation we label that with a T. The last letter in the names of models with radial flow is related to the speed of the radial 
flow. 



Models 


Infall type 




Threshold 


Radial inflow 


SII 


1 infall 


1.033 R (kpc) -1.27 Gyr 


no 


/ 


SIN 


1 infall 


4Gyr 


no 


/ 


S1IT4 


1 infall 


1.033 R (kpc) -1.27 Gyr 


4 Mopc-2 


/ 


S1IT7 


1 infall 


1.033 R (kpc) -1.27 Gyr 


7 Mopc-2 


/ 


S2IT 


2 infall 


1.033 R (kpc) -1.27 Gyr 


7 Mopc"^ (thin disk phase) 
4 Mopc"^ (halo thick disk phase) 


/ 


Rlll 


1 infall 


1.033 R (kpc) -1.27 Gyr 


no 


- 1 km s-i 


RlNl 


1 infall 


4 Gyr 


no 


- 1 km s-i 


RIIL 


1 infall 


1.033 R (kpc) -1.27 Gyr 


no 


linear inflow pattern 


RlITOl 


1 infall 


1.033 R (kpc) -1.27 Gyr 


4 Mopc-2 


-O.lkms-i 


R1IT05 


1 infall 


1.033 R (kpc) -1.27 Gyr 


4 Mopc-2 


- 0.5 km s-i 


R2ITV 


2 infall 


1.033 R (kpc) -1.27 Gyr 


7 Mopc~^ (thin disk phase) 
4 Mopc~^ (halo thick disk phase) 


variable inflow pattern 
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Fig. 4. Radial abundance gradient for oxygen. The black solid 
line refers to the one infall model without threshold and with 
the inside-out formation (SII). The red short dashed line repre- 
sents the model with a radial inflow of 1 km s"^ (RIH) and the 
blue long dashed line the best fit model using a variable veloc- 
ity for the radial inflow (RIIL). The data are the same of Fig. 

m 



old and inside-out is to produce a weak abundance gradient (~ 
-0.014 dex kpc-i) at variance with observations. In our model 
with a constant radial velocity of the flow, the abundance of 
oxygen at each galactocentric distance is lower than the val- 
ues found for the SIN model; this is due to the fact that the 
metals tend to be stored in the very central parts (R < 5 kpc) 
of the Galaxy and that we used a constant speed of the radial 
flow. In fact, Schonrich & Binney (2009) considered both the 
stellar and gas flows in a chemical evolution model consider- 
ing a model without inside-out formation. They were able to fit 




\ [kpc] 

Fig. 5. With the solid line we show the requested pattern of 
the velocity to reproduce the observed gradient as a function of 
the galactocentric distance for the RIIL model, with the dashed 
line for the R2ITV. In this plot we show the modulus of the ra- 
dial inflow velocity as a function of the galactocentric distance. 



the observed data including a variable speed of the gas flow by 
means of two free parameters, a situation where each ring may 
have its own velocity. 
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Fig. 6. Radial abundance gradient for oxygen. The black solid 
line refers to the one infall model without threshold and with 
the inside-out formation (SII). The red short dashed line rep- 
resents the model with a radial inflow of 1 km s"^ (Rlll). 
The data are mean values and relative errors from Cepheids 
by Andrievsky et al. (2002a,b). 
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Fig. 7. Radial abundance gradient for oxygen. The black solid 
line refers to the one infall model without threshold and 
with the inside-out formation (SII) with metallicity dependent 
yields, the green long dashed line is the case with solar yields . 
The red short dashed dotted line represents the model with a ra- 
dial inflow of 1 km s"^ (Rlll) with metallicity dependent yields 
and the blue long dashed line with solar yields. The data are 
mean values and relative errors from Cepheids by Andrievsky 
et al. (2002a,b). 



Fig. 8. Radial abundance gradient for iron. The black solid line 
refers to the one infall model without threshold and with the 
inside-out formation (SII). The red short dashed line represents 
the model with a radial inflow of 1 km s"^ (Rlll,) and the blue 
long dashed line the best fit model using a variable velocity for 
the radial inflow (RIIL). The mean values and relative errors 
from Cepheids by Andrievsky et al. (2002a,b) are reported with 
filled blue squares. 
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Fig. 9. Radial abundance gradient for oxygen. The black solid 
line refers to the one infall model without threshold and with 
the inside-out formation (SII). The red short dashed line rep- 
resents the model with a radial inflow of 1 km s"^ (RIH) us- 
ing shells 2 kpc wide and the blue long dashed line label the 
same physical parameters but running a model with shells 1 
kpc wide. The data are mean values and relative errors from 
Cepheids by Andriesky et al. (2005 a,b). 



We remind that also the SII model with inside-out forma- 
tion using the law described in eq. (O and no threshold is not 
able to reproduce the set of data for the outer part of the disk. 
Then we tried to find a best fit to the observed data varying the 
velocity of the gas flow. In Fig. [4] we report the Rlll model re- 
sults with a constant radial flow along the Galactic disk fixed at 
-1 km s"^ and model with an inflow velocity variable in space 
which results to be RIIL. We see that the model with the con- 
stant radial flow has an steeper gradient compared to SII but 
there are still problems in reproducing the outer region of the 
Galaxy. To reproduce the data we need a variable velocity for 
the inflow of gas. With the model RIITL we are able to fit quite 
well the data of HII regions and PNe assuming a linearly in- 
creasing flow velocity towards the outer regions as shown in 
Fig. (O. We recall that the case of a linearly increasing flow ve- 
locity was also studied in the work of Goetz & Koppen (1992). 



In Fig. [5] we show the pattern of the linear velocity we 
adopted to reproduce the observed gradient as a function of the 
galactocentric distance for the one-infall model. The range of 
velocities in modulus span the range 0-4 km s"^ in accordance 
with the results of Schonrich & Binney (2009). 

In Fig. Owe compare the models SII and Rlll with the data 
from Cepheids. We note that, since the observed gradient for 
these objects is almost flat, both models can easily reproduce 
it. 

We recall that in our calculations we consider the nucle- 
osynthesis prescription used in Fran9ois et al. (2004), therefore 
for oxygen we use metallicity dependent yields. To ascertain 
how the choice of metal dependent yields afl'ects radial flows, 
in Fig. [7] we run the same models of Fig. [6] but in the case of 
solar yields for oxygen. We note that the reference models as 
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Fig. 10. The surface gas density profile as a function of the 
radial distance for the reference one-infall model SlI drawn 
with the black solid line and for the best fit model RIIL. The 
green filled triangles represents the data of Dame et al (1993). 
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Fig. 11. Radial abundance gradient for oxygen. The black solid 
line refers to the one-infall model without threshold and with 
the inside-out formation (SlI). The red short dashed line rep- 
resents the model with a threshold of 4 M^pc'^ (S1IT4), the 
blue long dashed line the model with a threshold of 4 M^pc'^ 
combined to a radial inflow of 0.1 km s"^ (RlTOl), the green 
dotted line the model with a threshold of 4 MqPc~^ combined 
inflow of 0.5 km/s (R1T05). The data are the same of Fig.[2l 

well as the ones with radial flows, give exactly the same results 
for metal dependent and solar yields concerning the abundance 
gradients. This result is expected since the present time abun- 
dances in the ISM are the result of the global production of a 
given element; metal dependent yields vary with the metallicity 
but the global production of an element is the same as for solar 
metallicity yields since in both cases the solar abundances(i.e. 
the ISM abundances 4.5 Gyr ago) are reproduced. 

In Fig. [8] we show the results of models SlI, Rlll, and RIIL 
for iron. We note that iron shows a steeper gradient when com- 
pared with oxygen one, because of the difl'erent time-scales in- 
volved in the production of these two elements. 




Fig. 12. The SF rate normalized to the solar values as a function 
of the galactocentric distance for the one-infall model with no 
inside-out formation. The data are taken by Rana et al. (1991) 
and are reported with the magenta points and relative error b ars, 
with the red solid line we report the reference static model SIN 
whereas the dashed blue line is the model with a radial inflow 
of-1 kms"^ (RlNl). 




Fig. 13. The SF rate normalized to the solar values as a func- 
tion of the galactocentric distance for the one-infall model with 
inside-out formation. The data are taken by Rana et al. (1991) 
and are reported with the magenta points and relative error b ars, 
with the black solid line we report the reference static model 
SlI whereas the red short dashed line is the model with a radial 
inflow of -1 km s"^ (RIH) and the blue long dashed line is the 
best fit radial inflow model (RIIL). 
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Fig. 14. Radial abundance gradient for oxygen. The black long 
dashed line refers to the two-infall model S2IT. The blue solid 
line represents the the best 2 infall fit model using a variable 
velocity for the radial inflow (R2ITV). The data are the same 
of Fig. [2] 



Fig. 15. Radial abundance gradient for oxygen for the models 
SIN (black solid line), SlI (red dotted dashed line), and S2IT 
(green dotted line), and SI1T4 (magenta long dashed line). The 
data collection from Cepheids are reported. 



In our calculations we assume that: i) there is no inflow 
from the outer parts of the disk (R > 20 kpc), ii) our shells are 
2 kpc wide. 

The point ii) is consistent with the fact that on smaller 
scales other processes must be important. In fact, in Spitoni 
et al. (2008) we proved that galactic fountains can afl'ect abun- 
dance gradients only on scales smaller than 1 kpc. Here we test 
how the choice of a finer resolution aff'ects our results. In Fig. [9] 
the model Rlll is considered both for shells 1 and 2 kpc wide. 
We note that diff'erences arise only in the outer part of the disk, 
because of the adopted method for the implementation of radial 
inflows. It is worth noting that choosing a smaller shell does not 
afl'ect our results expect for distances R> 17 kpc and not in a 
substantial way as shown in Fig.[9j 

We pass now to analyze how the profile of the gas along 
the Galactic disk is modified with the variable velocity pattern 
used for the best fit of the abundance gradient. In the Fig. [TO] 
we compare the results of model SlI and RIIL for what con- 
cerns the gas surface density profile as a function of the radial 
distance. 

We note that both models underestimates the gas for dis- 
tances > 10 kpc and do not reproduce the peak and the decrease 
of gas for R < 5 kpc. This is probably due to the bar which is 
not considered in our model. 

We want also to stress that a model without a threshold in 
the SF has problems to reproduce the values in outer part of the 
disk. We also test a model in which we combine the eff'ect of 
a threshold in the SF plus a radial inflow. In Fig. [TT]is shown 
the "static " model S1IT4 with the "mixed models" considering 
both the threshold and the radial inflow. We fixed a threshold 
value at 4 M© pc"^ considering 2 constant speed radial inflows: 
the model RIlTOl with a velocity of 0.1 km s"^ and the model 
RI1T05 with a velocity 0.5 km s"\ respectively. A velocity of 
0.5 km/s is high enough to put system under threshold already 
at the radial distance of 16 kpc, whereas the model without 
radial flow but with a threshold of 4 M© pc"^ has no chemical 
evolution for R > 16 kpc. 

We want also to test how the radial inflows of gas afl'ect the 
SF rate along the disk of the Galaxy. In Fig. [12] we report the 



results for the SIN and the RlNl models, whereas in Fig. [13] 
the results with inside-out formation: SlI, Rlll and RIIL mod- 
els. In these Figs, we show the SF rate normalized to the solar 
value as a function of the galactocentric distance. Because of 
the uncertainties in this data set, as evident from the big error 
bars, we cannot draw firm conclusions. Having said that, we 
can however draw some considerations. In Fig.[T2]both models 
well fit the data for the SF rate, whereas in Fig. [T3]models SIl 
and Rlll underestimate the SFR/SFR© for R < 8 kpc. On the 
other hand, the model RIIL fits perfectly the SFR/SFR© ratio 
in the range 6-18 kpc. For galactocentric distances smaller 
than 6 kpc the RIIL model overestimates the mean SFR/SFR©, 
but our results are inside the 1 sigma error bar, because of the 
uncertainty in the data. We do not show our model results for 
galactocentric distances smaller than 4 kpc because in this re- 
gion the eff'ects of the central bar might be important and our 
model is built only for the disk and not for the bulge area. 

In this part we study the eff'ects of a radial inflow on the 
two-infall model described in Sect. 2. In Fig. [14] we see that our 
two-infall reference model S2IT is able to give rise to a steep 
gradient for oxygen but not enough to reproduce the data set. 
We also show the model R2ITV with a varying velocity pat- 
tern for the radial inflow velocity. In the case of the two-infall 
model we are not able anymore to fit the abundance gradient for 
oxygen using a simple constant relation between radial infall 
velocity and galactocentric distance, then we consider a pattern 
where the velocity starting from 17 kpc decreases monotoni- 
cally toward the center part of the Galaxy (see Fig. |5]l. 

From Fig. [5] we note that the range of velocities required 
to fit the data is quite similar to the one used for the one-infall 
model. 

In Fig. [15] we report the data from Cepheids compared with 
standard models without radial flows (SIN, SlI, SI1T4, S2IT). 
We see that all the models with inside-out formation show a 
rather shallow slope for the oxygen abundance gradient com- 
patible with the Cepheid data. In Fig. [16] we compare standard 
models without radial flows with all data sets considered in this 
paper. In conclusion, if we use only the data of the Cepheids, 
the abundance gradient is not so steep and it can be reproduced 
without any radial inflow. 
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Fig. 16. Radial abundance gradient for oxygen for the models 
SIN (black solid line), SlI (red dotted dashed line), and S2IT 
(green dotted line), and SI1T4 (magenta long dashed line). The 
data collection from Cepheids (blue squares ) and from HII 
regions (black circles) are reported. 



7. Conclusions 

In this paper we have studied the effects of radial inflows of gas 
on a detailed chemical evolution model. We also reported some 
results of the metals gradients obtained with "static" models 
considering different prescriptions for the infall law and star 
formation. Our main conclusions can be summarized as fol- 
lows: 

- If we consider the one-infall model with an inside-out for- 
mation (td varying with the radius) the obtained gradient 
without threshold is too fiat and the observational data are 
not reproduced especially in the outer part of the Galaxy 
disk. The slope for intermediate galactocentric distances 
can be reproduced if we consider the model with a thresh- 
old. However, since the surface gas mass density in the 
outer part of the disk is too small, there is no star forma- 
tion in this region and then no metals production and no 
chemical evolution. 

- Taking into account a constant td along the Galaxy disk 
equal to 4 Gyr, we find a fiat gradient in contrast with the 
work of Portinari & Chiosi (2000) where a gradient along 
the Galactic disk of 0.03 dex kpc"^ was found. 

- A constant radial inflow with speed of 1 km s"^ (not ex- 
pected theoretically) applied to the one-infall model with 
no threshold and no inside-out produces a weak abundance 
gradient ( — 0.014 dex kpc"^) and the metals tend to be 
stored in the very central parts of the Galaxy. 

- The required pattern of the the radial flow velocity varies 
linearly with the galactocentric distance and spans a range 
between and 4 km/sec. This conclusion holds for the one- 
infall model without threshold and with inside-out forma- 
tion. 

- Our "static" two-infall model with a threshold in the star 
formation and the inside-out formation is able to give rise 
to a steep gradient but not enough to reproduce the data set. 
The required pattern of the radial inflow velocity (in mod- 
ulus) to reproduce the observed gradient is quite similar to 
the one of the one-infall model. 



- If we use the data of the Cepheids the observed abundance 
gradient is not so steep and can be reproduced with our 
"static" reference models without any radial inflow. 

- As shown first in Gotz & Koppen (1992), then in Portinari 
& Chiosi (1999) and Colavitti et al. (2008), a variable effi- 
ciency in the star formation rate if coupled with an inside- 
out formation can also reproduce the abundance gradient. 
Since this variable efficiency in the Milky Way is observa- 
tionally motivated (Marcon-Uchida et al. 2010), it should 
not be disregarded as source of the abundance gradients. 
Moreover, as we have shown, to reproduce the abundance 
gradients by means of radial flows we need a specific vari- 
able speed of the flow, and at this moment there are not 
observational constraints for the speed of this flow of gas 
along the disk of the Galaxy. 

- Finally, we conclude that radial gas flows can be in prin- 
ciple important to reproduce the gradients along the disk, 
although an inside-out formation coupled with variable ef- 
ficiency of star formation and threshold in gas density can 
also well reproduce the data without radial flows. Probably, 
all these processes are present in the Galactic disk: in order 
to decide which of these processes would prevail in the for- 
mation of the disk, we need more detailed data both on the 
abundance, gas and star formation rate gradients, as well as 
data on high redshift disks. 

Acknowledgements. We thank the referee for the enlightening sug- 
gestions. We also thank G. Cescutti, and S. Recchi for many useful 
discussions. We aknowledge financial support from PRIN 2007 MUR 
Prot. No. 2007JJC53X-001. 

References 

Andrievsky, S.M., Bersier, D., Kovtyukh, V. V. et al 2002a, 
A&A, 384, 140 

Andrievsky, S.M., Kovtyukh, V.V., Luck R.E. et al 2002b, 
A&A, 381,32 

Asplund, M., Grevesse, N., Sauval, A. J., Scott, P., 2009, 
ARA&A, 47,481 

Berlin G., Lin C.C., 1996, Spiral structure in galaxies - A den- 
sity wave theory, MIT Press, Cambridge, Massachusetts 

Cescutti, G., Matteucci, R, Francois, P., & Chiappini, C. 2007, 
A&A, 462, 943 

Cescutti, 2008, A&A, 481, 691 

Chiappini C, Matteucci R, & Gratton R., 1997, ApJ, 477, 765 
Chiappini C, Matteucci R, & Meynet G., 2003, ApJ, 410, 257 
Chiappini C, Matteucci R, & Romano D., 2001, ApJ, 554, 
1044 

Colavitti. E., Cescutti, G., Matteucci R, Murante G., 2009, 
A&A, 496, 429 

Costa, R. D. D., Uchida, M. M. M., & Maciel, W. J. 2004, 

A&A, 423, 199 
Dame, T. M., 1992, AIPCP, 278, 267 

Deharveng, L., Peiia, M., Caplan, J., & Costero, 2000, 

MNRAS, 311,329 
Esteban, C, Garca-Rojas, J., Peimbert, M., Peimbert, A., Ruiz, 

M. T., Rodrguez, M., Carigi, L., 2005, ApJ, 618, 95 
Frangois P., Matteucci R, Cayrel R., Spite, M., Spite, R, & 

Chiappini, C, 2004, A&A, 421, 613 



Spitoni et al.: Effects of the radial flows 



11 



Gotz, M., Koppen, J., 1992, A&A, 262, 455 

Greggio L.,& Renzini A., 1983, A&A, 118, 217 

Kennicutt, R. C, Jr., 1989, ApJ, 344, 685 

Kennicutt, R. C, Jr., 1998, ApJ, 498, 541 

Iwamoto, K., Brachwitz, R, Nomoto, K., et al. 1999, ApJ 

Suppl. Ser. 125, 439 
Lacey CO., Fall M., 1985, ApJ, 290, 154 
Maeder, A., & Meynet, G., 1989, A&A, 210, 155 
Mannucci, F ., Delia Valle, M., Panagia, N., Cappellaro, E., 

Cresci, G., Maiolino, R., Petrosian, A., & Turatto, M, 2005, 

A&A, 433, 807 

Marcon-Uchida, M.M., Matteucci, F., Costa, R.D.D., 2010, 

WXiv: 1004.4139 
Martin, C. L., & Kennicutt, R. C., Jr., 2001, ApJ, 555, 301 
Matteucci, F., 2001, The Chemical Evolution Of The Galaxy, 

Kluwer Academic Publishers. 
Matteucci, F., & Fran9ois, R, 1989, MNRAS, 239, 885 
Matteucci M.F., Greggio L., 1986, A&A, 154, 279 
Mayor M., Vigroux L., 1981, A&A 98, 1 
Pagel, B.E.J.,& Tautvaisiene, G., 1995, MNRAS, 276, 505 
Portinari L., Chiosi C, 1999, A&A 350, 827 
Portinari, L., Chiosi, C, 2000, A&A, 355, 929 
Rana, N. C. 1991, ARA&A, 29, 129 

Rudolph, A. L., Fich, M., Bell, G.R., Norsen, T., Simpson, J.R, 
Haas, M.R., & Erickson, E.F., 2006, ApJS, 162, 346 

Scalo, J. M., 1986, FCPh, 11, 1 

Schmidt, M., 1959, ApJ, 129, 243 

Schonrich, R., & Binney, J, 2009, MNRAS, 396, 203 

Sommer-Larsen J., Yoshii Y., 1990, MNRAS, 243, 468 

Spitoni, E., Recchi, S., Matteucci, F., 2008. A&A, 484. 743 

Talbot, R., J., Jr., Arnett, W., D., 1973, 186, 51 

Thon R., Meusinger H., 1998, A&A, 338, 413 

Timmes, F.X., Woosley, S.E., Weaver, T.A., 1995, ApJS, 98, 
617 

van den Hoek L.B., & Groenewegen M.A.T. 1997, A&A 

Suppl., 123, 305 
Woosley S.E., & Weaver, T.A., 1995, ApJ, 101, 181 



